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Neutron capture by neutron-rich Li— and Be-isotopes plays a role in big- 
bang nucleosynthesis, especially in its inhomogeneous version and in the a— 
process occurring in supcrnovae. New reaction rates for 7,8 Li(n,7) 8,9 Li and 
9,10,11 Be(n,7) 10, n ' 12 Be have been consistently calculated using direct capture for 
the nonresonant part and the Breit-Wigner formula for the resonant part. The 
spectroscopic factors, spin/parity assignments and excitation energies of the final 
bound and initial resonant states have been taken from existing experimental data 
whenever possible. For unstable nuclei where this information is not experimen- 
tally available the shell model was used to determine these quantities. 



1 Introduction 

Evolution of the physical conditions of the universe, galaxies and stars can be 
described in terms of the increase or decrease of hundreds of elemental abun- 
dances of atomic nuclides. They originate from the primordial nucleosynthesis 
about fifteen billion years ago and the subsequent production/destruction cy- 
cle in stars and ejection into the intergalactic space. It is therefore inevitable 
and even fundamental to study the nuclear processes in several astrophysical 
sites for a deep understanding of the evolution of the universe. 

Cosmologically, the primordial nucleosynthesis provides a unique method 
to determine the average universal mass-density parameter f2 b ■ Although the 
homogeneous big-bang model for primordial nucleosynthesis predicts fls h% Q ~ 
0.04, X-ray observations of dense clusters have indicated that £Ib could be as 
large as < 0.15. Recent MACHO detections also suggest that there exist more 
baryons in our Galaxy than ever expected. There is clearly a serious poten- 
tial conflict between these observations and the theoretical prediction in the 
homogeneous big-bang model. The situation is even crucial if high deuterium 
abundances, which were detected in Lyman-a absorption systems along the 
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line of sight to high red-shifted quasars, are presumed to be primordial. On the 
other hand, an inhomogeneous big-bang model, which allows inhomogeneous 
baryon density distribution, can predict ft|o ~ 0.1 — 0.2. Among the pos- 
sible observable signatures of baryon inhomogeneous cosmologies are the high 
abundances of heavier elements than lithium such as beryllium and borona. 
In an environment of baryon inhomogeneous distribution, neutrons can easily 
diffuse out of the fluctuations to form high density proton-rich and low density 
neutron-rich regions, where a lot of proton/neutron-rich radioactive isotopes 
can help produce the intermediate-to-heavy mass elements. 

Another astrophysical site where the neutron-rich isotopes may play a sig- 
nificant role in nucleosynthesis is the a-process occurring in supernovae. The 
nucleosynthesis in the high-entropy bubble is thought to proceed as follows. 
Due to the high temperature, the previously produced nuclei up to iron will 
be destroyed again by photodisintegration. At temperatures of about 10 10 K 
the nuclei would be dismantled into their constituents, protons and neutrons. 
At slightly lower temperatures one is still left with a-particles. During the 
subsequent cooling of the plasma the nucleons will recombine again, first to 
a-particles, then to heavier nuclei. Depending on the exact temperatures, den- 
sities and the neutron excess, quite different abundance distributions can be 
produced in this a-rich freeze-out (sometimes also called a-process) . Temper- 
ature and density are dropping quickly in the adiabatically expanding high- 
entropy bubble. This will hinder the recombination of a-particles into heavy 
nuclei, leading in some scenarios to a high neutron density for an r-process, at 
the end of the a-process after freeze-out of charged particle reactions. 

These astrophysical motivations have led us to critically study the the role 
of radiative neutron capture reactions by neutron-rich Li- and Be-isotopes 
theoretically in explosive nucleosynthesis. Since it is the focus in recent years 
to study the the nuclear reactions dynamics by the use of radioactive nuclear 
beams, our theoretical studies are also being tested experimentally. 

2 Calculation of Radiative Capture Cross Sections 

Nuclear burning in explosive astrophysical environments produces unstable 
nuclei which can again be targets for subsequent reactions. In addition, it 
involves a very large number of stable nuclei which are not yet fully explored 
by experiments. Thus, it is necessary to be able to predict reaction cross 
sections and thermonuclear rates with the aid of theoretical models. 

In astrophysically relevant nuclear reactions two important reaction mech- 
anisms take place. These two mechanisms are compound-nucleus reactions 
(CN) and direct reactions (DI). The reaction mechanism and therefore also 
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the reaction model depends on the number of levels in the CN. If one is con- 
sidering only a few CN resonances the R-matrix theory is appropriate. In 
the case of a single resonance the R-matrix theory reduces to the simple phe- 
nomenological Breit-Wigner formula. If the level density of the CN is so high 
that there are many overlapping resonances, the CN mechanism will dominate 
and the statistical HF-model can be applied. Finally, if there are no CN reso- 
nances in a certain energy interval the DI mechanism dominates and one can 
use DI models, like Direct Capture (DC). 

In the case of a single isolated resonance the resonant part of the cross 
section is given by the well-known Breit-Wigner formulantl: 

. . _ nh 2 (2 j + 1) r in r ou t . . 

ai[ ' ~ 2fiE (2j p + 1) (2j t + 1) (£ r _ E f + Eik ' U 

where J is the angular momentum quantum number and E T the resonance 
energy. The partial widths of the entrance and exit channels are Ti n and r out , 
respectively. The total width r to t is the sum over the partial widths of all 
channels. One important aspect is that the particle width T p can be related 
to spectroscopic factors S and the single-particle width r s . p . bjoEl 

r p = c 2 sr s . p . , (2) 

where C is the isospin Clebsch-Gordan coefficient. The single-particle width 
r s p . can be calculated from the scattering phase shifts of a scattering potential 
with the potential depth determined by matching the resonance energy. 

The-jaonresonant part of the cross section can be obtained using the DC 
modeBOB: 

= C 2 S c o? c . (3) 

c 

The sum extends over all bound states in the final nuclei. The DC cross 
sections cr? are essentially determined by the overlap of the scattering wave 
function in the entrance channel, the bound-state wave function in the exit 
channel and the multipole transition-operator. 

The total cross section can be calculated by summing over the resonant 
(Eq. |l|) and nonresonant parts (Eq. |J) of the cross section (if the widths of 
the resonances are broad, also an interference term has to be added). For 
both parts the spectroscopic factors have to be known. They can be obtained 
from other reactions, e.g., the spectroscopic factors necessary for calculating 
A(n,7)B can be extracted from the reaction A(d,p)B. The 7-widths can be ex- 
tracted from reduced electromagnetic transition strengths. For unstable nuclei 
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where only limited or even no experimental information is available, the spec- 
troscopic factors and electromagnetic transition strengths can also be extracted 
from nuclear structure models like the shell model (SM). 

The most important ingredients in the potential models are the wave func- 
tions for the scattering and bound states in the entrance and exit channels. 
This is the case for the DC cross sections <7^? c in Eq. || as well as for the cal- 
culation of the single-particle width in Eq. g. For the calculation of these 
wave functions we use real folding potentials which are given byurQ 

V{R) = \V F {R) = \ J J p a {r 1 )p A {v 2 )v cS {E, Pa ,p A ,s)&r 1 &r 2 , (4) 

with A being a potential strength parameter close to unity, and s = |R+r2 — ri | , 
where R is the separation of the centers of mass of the projectile and the 
target-nucleus. The density can been derived from measured charge distribu- 
tionsE3 or from nuclear structure models (e.g., Hartree-Fock calculations) and 
the effective nudeon-nucleon interaction v D g has been taken in the DDM3Y 
parametrizatioriJ. The imaginary part of the potential is very small because of 
the small flux into other reaction channels and can be neglected in most cases 
involving neutron capture by neutron-rich target nuclei. 



3 Reaction Rates for Li and Be Isotopes 

The parameters for the resonant and nonresonant contributions to the reaction 
rates are listed in Tables [j] and ||, respectively. In the tables we give experimen- 
tal values if available. Otherwise the excitation energies, spectroscopic factors, 
neutron- pad 7-widths were calculated with the shell model. We used the code 
OXBASHliJ for the calculations. For normal-parity states we employed the in- 
teraction (8-16)POT of Cohen and Kurathli3. For noiinormal parity states we 
used the WBN interaction of Warburton and Browr£3. 

With Eq. [I] the resonant reaction rate can be derived as 

iV A (< 7V ) i = 1.54 x l0 5 p- 3/2 T~ 3/2 (5) 
^7* exp(-11.605£ r /T 9 ) cm 3 mole" 1 s" 1 , 

i 

where Tg is the temperature in 10 9 K, E r the resonance energy in the cm. sys- 
tem (in MeV), and the resonance strength 0^7 (in eV) is given by 

2J+ 1 

U1 = (2j p + i)(2 Jt + i) r tot ' [b) 
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Table 1: Resonance parameters 



Reaction 


E x 


E n 




r„ 


r 7 


W7 




(MeV) 


(MeV) 




(eV) 


(eV) 


(eV) 


Y Li(n, 7 )«Li 


2.26 


0.227 


3+ 


3.1 x 10 4 


0.07 


0.061 


8 Li(n, 7 ) 9 Li 


4.31 


0.247 


5/2- 


1 x 10 5 


0.11 


0.066 


9 Be(n,7) 10 Be 


7.371 


0.559 


3- 


1.57 x 10 4 


0.661 


0.578 




7.542 


0.73 


2+ 


6.3 x 10 3 


0.814 


0.509 



The partial widths of the entrance and exit channel, T m and r ou t, are in the 
case of (n,7)-reactions the neutron- and 7-widths. Since the neutron width 
is usually much larger than the 7-width, the total width r to t is practically 
identical with the neutron-width. 

In Table [l]we list the excitation energies, resonance energies, neutron- and 
7-widths and the resonance strengths of the resonances. 

The nonresonant capture cross section is parametrized as 

O'nr(E) — A/s/E + Bs/E — CE D , (7) 

with [A] = fib MeV 1 / 2 , [B] = ^bMeV" 1 / 2 , and [C] = fib MeV " D . The param- 
eters A,B,C and D are listed in Table ||. Using this equation, we obtain for 
the reaction rate 

(<™) nr = ^836.565,4/!- 1/2 + 108. 130 Bfr 1/2 T g 

-277.097^- 1 / 2 ^^T ! f +1/2 )cm 3 s- 1 mole- 1 , (8) 

where fi ist the reduced mass in units of the atomic mass unit and V{z) is the 
Euler gamma function. 

The total reaction rate is given as the sum of the resonant (Eq. || and 
nonresonant (Eq. |^) part. 

3.1 7 Li(n,j) 8 Li 

The cross section of the reaction 7 Li(n,7) 8 Li is well known (see, e.g.,0). The 
cross section is dominated by s-wave capture to the 8 Li ground state and 
a resonance at 227_keV neutron energy. Using the spectroscopic factors of 
Cohen and KurathEJ yields a thermal cross section of 8.2 x 10 -2 b, which is a 
factor 1.8 higher than the experimental value of 4.54 x 10 -2 b. The shell model 
calculation is purely p-shell and does not include excitations to other oscillator 
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Table 2: Parametrization of the nonresonant cross section (see text). 





A 


B 


c 


D 


Cnr(Mb) at 30keV 












This 


Rauscher 












work 


et aM 


7 Li(n, 7 ) B Li 


6.755 a 








39.000 




8 Li(n, 7 ) 9 Li 


2.909 








16.795 


30.392 


9 Be(n,7) 10 Be 


1.147 a 


11.000 


6.815 


0.962 


8.294 


6.622 


10 Be(n,7) n Be 


0.132 


24.000 


15.725 


0.914 


4.281 


3.943 


n Bc(n,7) 12 Bc 




7.000 


4.851 


0.887 


0.996 


2.373 



a extracted from experimental thermal cross sectioi 



shells. Therefore the spectroscopic amplitude of 0.977 for a p 3 / 2 -transition to 
the ground state of 8 Li might be too high. 

For the resonance, however, we find excellent agreement between calcu- 
lation and experiment. The calculated width — using the folding potential 
and spectroscopic amplitudes from Cohen and Kurathlia — is 28.9 keV, almost 
identical to the known value of 31 ± 7keV. 

3.2 8 Li(n,j) 9 Li 

The resonance at 247 keV is a 5/2" stateB. With a total width of 100 keV 
the resonance strengthjJs determined by the 7-width which was previously 
estimated with 0.56eVlia. A shell model calculation yielded a width T 7 = 
0.11 eV. Therefore the resonance strength is a factor 5 smaller than previously 
assumecO. 

The calculated thermal cross section, resulting from s-wave capture to the 
ground state and first excited state in 9 Li, is 1.94 ^f-\0 2 b and is smaller than 
the value of 3.51 x 10~ 2 b given by Rauscher et al£3. 

3.3 9 Be(n,j) w Be 

Likc_m the reaction 7 Li(n,7) 8 Li the spectroscopic factors of Cohen and Ku- 
ratbli3 are a little too high. The thermal cross section is dominated by the 
transition to the 10 Be ground state with a theoretical spectroscopic factor of 
2.36. With this value the calculated thermal cross section is 1.06 x 10~ 2 b, com- 
pared to experimental cross section of 7.6 x 10 -3 b. With the spectroscopic 
factor given by MughabghahcJ of 1.45 the calculated cross section would be 
close to the experimental value. For high temperatures the p-wave capture to 
excited states has to be taken into account. 
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Two resonances are known at 559 keV and 730 keV. The total widths are 
known experimentally. We have calculated the 7-widths which were only es- 
timated previously. Both resonance strengths are larger than the previous 
estimates, for the 559 keV resonance the enhancement is one order of magni- 
tude. 

With the higher resonance strengths and .the p-wave contribution the re- 
action rate is clearly higher compared to Ref. t3. 

3.4 w Be(n, 1 ) 11 Be 

Cross section and reaction rate of this reaction were recently-determined exper- 
imentally with the help of the inverse Coulomb dissociatior£H. They supported 
their experimental values by a direct capture calculation. In order to repro- 
duce the experimental data they enhanced the spectroscopic factors to the 
11 Be ground state by 20%. Our calculation pCpnfirms the results. Using the 
spectroscopic factors from the (d,p)-reactionE3 the calculated cross section is 
a little smaller than the experimental. The results are grossly different from 
the rate given in Ref. Ej. 

3.5 ll Be(n,7) 12 Be 

There is no resonant contribution to the reaction rate. The transition is a p— 
wave capture from the l/2 + ground state of 11 Be to the ground state of 12 Be 
and the + state at 2.7 MeV excitation energy, while the transition to the 2 + 
state at 2.1 MeV is negligible. 

4 Discussion 

The new reaction rates could change the reaction flow in the inhomogeneous 
big bang nucleosynthesis. The smaller rate for 8 Li(n,7) 9 Li could mean that the 
main reaction flow will proceed through the reaction 8 Li(a,n) 11 B. The higher 
rate for 9 Be(n,7) 10 Be might give more importance to this reaction. Detailed 
network calculations with the new rates are planned for the near future. 
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